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ABSTRACT 

We find a bright (Ltjv = 2.5L* =6 ) Lyman alpha emitter (LAE) at redshift z — 
5.719 at a projected distance of 79 physical kpc from a strong triply ionized carbon 
(C IV) absorption system at redshift z = 5.7238 previously reported in the spectrum of 
the Zcm = 6.309 QSO SDSS J1030+0524. This is the highest redshift galaxy-absorber 
pair detected to-date, supporting the idea that galaxy-wide outflows were already in 
place at the end of the epoch of reionization. The proximity of this object makes it 
the most likely source of metals, consistent with models of outflows at lower redshift 
where significant observational evidence relates metal absorption systems with galaxies 
hosting outflows. 

In a typical outflow scenario, a wind of 200 km s , active since the universe was 
only 0.6 Gyr old (z ~ 8.4), could eject metals out to 79 kpc at z = 5.719. Although 
the origin of metals in the intergalactic medium (IGM) is still under debate, our 
results are consistent with predictions from cosmological simulations which reproduce 
the evolution of the cosmic density of C iv from z ~ 6 to the present day based on 
outflow-driven enrichment of the IGM. 

We also report two more Lya emitters in this field, at z — 5.973 ± 0.002 and 
z = 5.676 ± 0.002 respectively, the former confirming the original identification by 
Stiavelli et al. Our results suggest that the colour cut typically used to identify i- 
dropouts (1775 — zg5o > 1.3) misses a non-negligible fraction of blue galaxies with faint 
UV continuum at z ^ 5.7. 

Key words: early universe, galaxies: high redshift, galaxies: intergalactic medium, 
galaxies: distances and redshifts. 



1 INTRODUCTION 

Medium and high resolution spectroscopy has revealed 
the presence of metals in the high redshift intergalactic 
medium (IGM) through the detection of absorption sys- 
tems in the spectra o f background quasi stellar objects 



(QSOs) (e.g. ISongailal l200ll 
20061: IBecker et al 



iPettini et al. I 120031: ISimcod 



20061: IBecker. Rauch &. Sargent! 120091: 
Rvan-Weber et all 120091 : IBecker et all l201ll : ISimcoe et all 
201 If ). However, the enrichment source of these absorption 
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systems is still poorly constrained. There is a debate as 
to whether these intergalactic metals originated in the re- 
cent past from nearby galaxies driving large-scale outflows 
of their interstellar media, or whether the IGM may have 
received a 'blanket' enrichment at much earlier times, when 
the Universe became reionised by the first stars and galaxies. 



IPorciani fc Madaul (|2005h favour the "pre-galactic en- 
richment" scenario, arguing that the cross-correlation func- 
tion between Lyman break galaxies (LBGs) and Civ ab- 
sorbers at z ~ 3, as well as the mean absorption line den- 
sity, can be explained by enriched supernova (SN) driven 
wind s from dwarf galaxies at z > 6 jMadau. Ferrara fc Reesl 
l200ll) . These LBG progenitors could generate metal bubbles 
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~ 100 kpc (comoving) in size within the same high-density 
regions of the IGM where LBGs would later form. 

On the other hand, the connection between star 
forming galaxies and metal absorption features is now 
firmly established at redshifts 1 ^ z £ 3, strongly support- 
ing the idea that galactic outflows, or 'superwinds' as 
they are often referred to, are one of the main sources of 
metal enrichment of the IGM, a t least at these redshifts 
(|Steidel et all |2010| ; 1 imcoel l201ll . and references therein) . 
Evidence for such outflows is provided by the nearly uni- 
versal blueshift of interstellar absorption lines observed 
against the nuclei of g alaxies undergoing vig o rous star for- 



mation episodes (e.g. Hcckma n et al 



200(1 iPettini et all 



l200ll ; IShaplev et al.ll2003l : ISteidel etaTl^OlOh . with the out- 
flowing gas reaching velocities of up to ~ 1000 km s -1 rela- 
tive to systemic (e.g. [Pettini et al.ll2002l ; lOuider et aT]|2009l ; 
iDessauees-Zavadskv et al.ll20ld . and reference therein). In a 
recent study of galactic outflows at z ~ 1.4, IWeiner et al.l 
(2009) found that wind velocities increase with mass and 
star formation rate (u w ind oc SFR ' 3 ). A similar trend 
(fwind oc SFR ' 35 ) had previosuly been shown to hold 
in ultra luminous inf rared ga l axies (ULIR Gs) at redshift 
z = 0.042 - 0.16 bvlMartinl <l2005h: both iMartinl l|2005l ) 
and iRupke. Veilleux fc Sanders! (|2005l ) found that 70-85% 
of ULIRGS at low redshifts [z < 0.5) drive such energetic 
outflows. 

The common occurrence of outflows in galaxies under- 
going bursts of star formation and the obvious importance 
of this phenomenon for galaxy evolution and the enrichment 
of the IGM are some of the reasons why superwinds are at- 
tracting a great deal of theoretical interest. One of the goals 
of current theoretical work is to provide a physical explana- 
tion for the somewhat puzzling redshift evolution of the co- 
moving density of Civ ions in the IGM, flQjy^z). Looking 
back in time, Q(j l y(z) appears to be approximately con- 
stant during a ~ 3Gyr interval between z ~ 1.5 and 5, 
and then drops by a factor of ~ 3.5 over a period of only 
~ 300 Myr from z = 4.7 to 5 . 8 dRvan- Weber et all 120091 ; 
iBecker. Rauch fc Sargentll2009l; Isimcoe et alj|201ll) 
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collaborators have explored different wind scenarios using 
cosmological hydrodynamic simulations which explicitly in- 
corporate outflows from starburst galaxies and concluded 
that momentum-driven winds give the best agreement with 
the observations. Interestingly, in their simulations the out- 
flows travel to similar distances -Rturn = 80 ± 20 physical 
kpc at all redshifts they considered. If this is the case in 
reality, it would imply that at early epochs the metal en- 
richment by star-forming galaxie s affects larger volum es of 
the Universe than at later times. iTescari et al.l (|2010Tl have 
also highlighted the importance of feedback from galactic 
winds for understanding the properties of metal absorption 
systems in the IGM. Among the models they considered, 
those without galactic feedback could not reproduce the ob- 
served evolution of the Qq iv (z), while momentum-driven 
wind simulations gave the best ma tch with observations, 
in agr eement wit h the proposal by Oppenheimer fc Pavel 



(2006). Similarly. ICen fc Chisarl < |201Cf ) concluded that the 
dominant mechanism responsible for the metal enrichment 
of the IGM is star formation feedback from the analysis of 



hydrodynamical simulations tuned to reproduce the star for- 
mation history of the universe from z — 6 to z — 0. In these 
simulations, C IV and O VI are mostly located in regions with 
temperature T ^ 2 x 10 4 K that have been swept-up by the 
shocks produced by galactic-scale winds; absorption systems 



with column densities N, 



Civ 



10 -10 1 



ated with galaxies of luminosities L L* . 

While it is generally agreed that galactic outflows are 
the most likely mechanism for enriching the IGM, the obser- 
vational evidence for such outflows has so far been limited 
mostly to galaxies at z 4. At higher redshifts, even L* 
galaxies are too faint to record individual spectral features 
- other than the Lya emission line - with sufficient precision 
to measure kinematic offsets betw een the stars and the out- 
flowing gas in individual galaxies. IVanzella et al.l (|2009T ) did 
find such offsets in stacked spectra of z = 5 and 6 galaxies 
selected from the Great Observatories Origins Deep Survey 
and indeed concluded that: "powerful, large-scale winds are 
common at high redshift". Nevertheless, it is of great in- 
terest to explore the outflow phenomenon in more detail at 
the highest redshifts, as it holds the key not only to the 
initial stages of the chemical evolution of the Universe, but 
also presumably to the escape of Lyman continuum photons 
responsible for reionizing the IGM at z > 6. 

If the general picture assembled from the above cos- 
mological simulations is correct, most Civ absorption sys- 
tems at z > 5 should be on their first outward tra- 
j ectory and should be found close t o the parent galaxy 
l|Oppenheimer. Dave fc Finlatorl 120091 ') . Then, a key obser- 
vational test of these models is to identify the star-forming 
galaxies close to known C IV absorption systems in the spec- 
tra of the highest redshifts QSOs. 

One of the strongest C IV absor bers in the high- 
z surv eys by iRvan- Weber et all ||2009h and ISimcoe et al.l 
|201ll) is at z abs = 5.7238 towards the z em = 6.309 
QSO SDSS J1030+0524, with rest frame equivalent 
widths W (A1548) = .65 A and IF (A1550) = 0.41 A 
jRvan- Weber et alj|2009h . The field of this QSO was im- 
aged in the 1775 and zsso bands with the Advanced Camera 
for Surveys (ACS) on the Hubble Space Telescope (HST) by 
IStiavelli et atl |2005h . These authors reported a statistical 
excess of galaxies photometrically classified as i-dropouts 
(i.e. with 2775 — zsso > 1.3), and hence pot entially located at 



z 5.5. This result was later confirmed bv lKim et al.l (|2009l ) 
as part of a larger study of i-dropouts around five QSOs at 
z ~ 6. In the latter work, 14±4 i-dropouts were found in the 
field of the QSO SDSS J1030+0524 (~ 11.3arcmin 2 ) com- 
pared to 8±3 expecte d from a random distr ibution of galax- 
ies at these redshifts |Giavalisco et al.1 l2004). In the absence 
of spectroscopic redshifts for all b ut one of the ii-dropo uts - 
a Lya emi t ter at z = 5.970, both lStiavelli et al l (|2005h and 
iKim et al] l|2009T l speculated that the excess galaxies may 
be part of the same overdensity that includes the QSO. 

However, it may also be the case that some of the candi- 
date galaxies are in fact associated with the foreground C IV 
absorber, particularly as two of the i-dropouts lie within 
~ 15.5arcsec of the QSO. This separation corresponds to 
a transverse distance of ~ 90 kpc (physical) from the QSO 
line-of-sight at z — 5.7; at 2 = 2-3 this is the impact pa- 
rameter within which strong C IV absorption i s found to be 
assoc iated with actively star-forming galaxies (Steid eTet al.l 
2010). Thus, this field is an ideal location for testing whether 
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galactic outflows are indeed responsible for producing C IV 
absorbers at z > 5.5, and how the properties of such out- 
flows compare with those at lower redshifts. 

Motivated by these considerations, we have performed 
deep spectroscopic observations of the two galaxies closest to 
the sight-line to the QSO SDSS J1030+0524, and of the ad- 
ditional galaxy alre ady known to be a Ly q emitter from the 
discovery paper bv lStiavelli et all (|2005T ). This paper is or- 
ganised as follows. In Section 2 we describe our spectroscopic 
observations as well as our reass essment of the photom etry 
of the ACS images obtained by IStiavelli et all (|2005l ). We 
present our results in Section 3 and discuss their implica- 
tions in Section 4. We summarise our conclusions in Section 
5. Throughout this work we use AB magnitudes and assume 
a fiat universe with Ho = 70 km s _1 Mpc -1 , fi m = 0.3 and 

n x = 0.7. 



into consideration that the ESO pipeline was developed 
specifically to reduce any systematic effect from the tele- 
scope and the instrument, we resumed the reduction pro- 
cess of ESO pipeline output files using iraf. Given that the 
galaxies are very faint, the sky subtraction residuals in each 
slit were removed with the task BACKGROUND, using a first 
order polynomial, before the extraction of individual spec- 
tra. Automatic cosmic ray cleaning was carried out at this 
stage, with a few cosmic-ray affected pixels requiring man- 
ual cleaning. Finally, after flux calibrating each spectrum, 
we combined the spectra of each galaxy in a weighted aver- 
age, where the weight is the median flux value of the 8-10 
pixels centred on the Lycv emission line, which is the only 
spectral feature detected. The results of our spectroscopic 
analysis are presented in Section 3. 



2.2 Photometry 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Spectroscopic Observations 

We obtained medium resolution spectra of three galaxies 
in the field of the QSO SPSS J1 030+0524 photometrically 
identified bv lStiavelli et all {2005) to be at at z > 5.5. Two of 
the galaxies were selected as being closest to the QSO sight- 
line, at transverse distances of ~ 90 kpc (assuming that they 
are at z = 5.5-6). The third galaxy, J103024.08+052420.4 
is the only object in the field for which lStiavelli et all (2005) 
were able to measure a spectroscopic redshift, z = 5.970; we 
re-observed this source for comparison purposes. 

The spectra were obtained with FORS2 on UT1 of the 
ESO VLT, operated in multi-object spectroscopic mode us- 
ing 0.7" slits with the GRIS_600z grism which covers the 
wavelength range from 7370 to 11000 A. The resolution is 
FWHM = 4.9A (or 172.5 km s _1 at 8500 A), sampled with 
three pixels. A total of 36 exposures of 900 s each were ac- 
quired between 2007 April 24 and 2007 June 17; two of these 
were of very low signal-to-noise ratio and were excluded from 
the final stack of combined spectra. The total integration 
time of the final image was 30 600 s. 

The observations were carried out in service mode and 
processed by ESO using automatic pipelines. The data prod- 
ucts consist of a sky-subtracted 2D image for each 900 s ex- 
posure, wavelength calibrated and corrected for bias, over- 
scan, dark current and flat field. Given the challenging na- 
ture of these faint object observations, we also processed the 
images independently of the ESO pipeline (starting from 
the sky subtraction stage), using conventional iraf tasks. 
Specifically, we combined the two-dimensional slit image of 
each object with imcombine and we found a single emission 
line centered along the spatial axis (slit direction) at the 
position where each object is expected to be. Then, we ex- 
tracted individual exposures with APALL and combined them 
in observational blocks of three exposures with SCOMBINE 
prior to the flux calibration. Once each block was flux cali- 
brated (using the provided spectroscopic standard star spec- 
tra and iraf tasks standard, sensfunc and calibrate), 
we averaged all the blocks and confirmed the presence of an 
emission line uncontaminated by sky residuals in the final 
combined spectrum of each object. 

Having performed this independent check, and taking 



The existence of an overdensity of galaxies at z > 5.5 in the 
fiel d of the z = 6. 3 09 QS O SDSS J1030+0524 first reported 
by^ 
by 



Stiavelli et all (120051 ) has more recently been confirmed 
Kim et al.l ( 20091 ). This was not unexpected, given that 



both studies used the same ACS/ HST images, the same 
colour section criterion and the same S/N limit. It is there- 
fore somewhat surprising that two out of the three galaxies 
we targeted, originally communicated to us by M. Stiavelli, 
[in lStiavelli et all l|2005l ) no information was provided on in- 
dividual galaxy positions, magnitu des and colours], a re not 
included in the later compilation bv lKim et all |2009). Hav- 
ing confirmed spectroscopically that both objects are indeed 
at z > 5.5 (at z = 5.676 and z — 5.719 respectively; see Sec- 
tion 3), we decided to remeasure ourselves the magnitudes 
^850 and colours 1775 — zgso of all three galaxies considered 
here. 

The data retrieved from the public HST archived] con- 
sist of seven images in the 1775 band with a total exposure 
time of 5840 s, and nine images in the zsso band with a total 
exposure time of 11300 s. We ran multidrizzle on the re- 
duced images (fit files) with a linear size of output pixels in 
arcseconds per pixel of pixscale = 0.6. The size of a pixel in 
units of the input pixel size was ptxfrac = 1.0 and inverse- 
variance maps were used as weighting image (finaLwht_type 
= ivm). Distortion corrections were computed with the most 
updated distortion coefficient tables available from the Mul- 
timission Archive at the Space Telescope (MAST). 

We used SEXTRACTOR version 2.5.0 to calculate the 
magnitudes and colours of the three galaxies in our sam- 
ple in order to have a consistent comparison with previ- 
ous works in the same field. Photometry was performed in 
pseudo-dual image mode in the zsso band and dual image 
mode in the 1775 band with the Z850 as detection image. The 
weight map from multidrizzle was used in the detection 
image and the corresponding exposure time map was used in 
the measurement images. We adopted MAG_AUTO as the 
total magnitude of each source. Zeropoint magnitudes were 
obtained from the Space Telescope Science Institute web- 
sitt[3 updated on May 19th, 2009. The values corresponding 
to F775W and F850LP filters for the date of observation 



1 http://archive.stsci.edu 
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are 25.67849 and 26.86663 respectively. Regard ing extinc- 
tion c orrection, the same values determined by iKim et aU 
(2009) were used: 0.048 and 0.036 for the 2775 and 2850 mag- 
nitudes respectively. Colours give in Table [1] were computed 
using isophotal photometry (MAGJSO). 

Owing to correlated noise introduced by multtdrizzle 
in the final image, the errors in magnitude and flux determi- 
nation provided by SEXTRACTOR are incorrect. To overcome 
this problem, the first step is to calculate the following cor- 
rection factor, A, valid in the case pixscale (s) ^ pixfrac (p) 



A = (-_ 



— f- 

p 3p 2 

Then, we applied the following correction to the error in the 
flux estimation and the error in the magnitude 
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where / crr is the flux error estimated by SEXTRACTOR, / is 
the measured flux and g is the gain. 

In Table \T\ we list the magnitudes and colours of the 
three galaxies corrected for correlated noise a nd galactic 
extinc tion. For the one galaxy in common with Kim et al.l 
l|2009l ). J103024.08+052420.41, we fin d z 850 = 25 94 ± 0.05 
and 1775 — 2850 = 2.25; for comparison IKim et al.l (|2009T ) re- 
ported 2850 = 25.74 and 1775 — 2 8 5o = 2.f2 (object A13 in 
their Table 3). 



3 RESULTS 

3.1 Lyman alpha detections 

As shown in Figs. 1, 2, and 3, we detect a single emission 
line in each of the three galaxies observed. We interpret these 
features as Lya emission lines at redshifts 2 em = 5.7-6.0, for 
the following reasons. 

As discussed in more detail below (Section 1 3 - 3 P - the 
(«775 — 2850 ) colours of the galaxies together with their faint 
2850 magnitudes are indicative of high redhifts, 2 > 5.4. 
The main contaminants in this colour regime are galaxies at 
2 = 1-2; in this case, the emission lines we have detected 
could conceivably be [0 11] AA3726, 3729. at z em = 1.2-1.3. 
However, the [O n] doublet, with a wavelength separation 
of 2.8 x (1 ± 2) ~ 6.2 A should be resolved in our spec- 
tra, which have a resolution of 4.9 A FWHM. Furthermore, 
even if the individual doublet lines were not fully resolved, 
we expect the width of the features to be FWHM ^ 8 A, 
broader than FWHM ~ 4.9 and 6.0 A we measure for the 
two newly discovered emission lines reported here, in Tar- 
gets 2 and 3 respectively. The emis sion line in Target 1 was 
independently identified as Lya bv lStiavelli et all (2005) Oil 
the basis of its asymmetric profile which is characteristic 
of resonantly scattered L ya emission in galax ies undergoing 
large-scale outflows (e.g. iQuider et aT]|2009l ). We conclude 
that the most likely interpretation is that all three features 
are Lya emission lines at redshifts 2 cm = 5.7-6.0, as detailed 
below. 

We measured the values of 2 om for each galaxy from 
the peak of the emission line; we estimate the uncertainty of 




03024.08 + 052420.4 
= 5.973 ± 0.002 
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Figure 1. FORS2 spectrum of J103024.08+052420.41. Top: 2D 
image of the spectrum. The x-axis is wavelength and the y-axis is 
the spatial direction along the slit. Bottom: ID spectrum, with the 
dotted line showing the lcr error spectrum. Grey regions indicate 
positions of all the sky lines present (with different intensities). 
In both 2D and ID spectra, an emission line is clearly detected 
at 8477 A; we identify this feature as Lya at z = 5.973. 



the redshifts to be ±0.002 from consideration of the seeing 
and spectral resolution of our data. However, it should be 
borne in mind that the Lya forest could absorb a significant 
fraction of the emission on the blue wing of the line, leading 
to an overestimate of the redshift. 

We also measured the line flux, from which we calculate 
L-Ly a , the Lya luminosity in our cosmology. This in turn 
can be used to estimate the star formation rate, SFRL yQ , 
using Kennicutt's (1998) calibration of SFRhcv and case B 
recombination 



SFR LyQ (M yr- 1 ) = 9.1 x 10- 43 L Lya (erg S - 



(1) 



The above conversion assumes a standard Salpeter 
(1955) stellar initial mass function (IMF); adopting the 
Chabrier (2003) IMF which has proportionally fewer low 
mass stars would reduce the values of SFRL ya by a factor of 
~ 1.8. On the other hand, it is well known that the uncor- 
rected Lya emission line luminosity can lead to significant 
underestimates of the SFR - not only because of overlapping 
absorption by the Lya forest, but especially as a consequence 
of resonant scattering. Multiply scattered Lya photons can 
be absorbed by even small quantities of dust and can diffuse 
spatially over extended halos greater than the projected size 
of the spectrograph slit (Steidel et al. 2011). We now briefly 
discuss each galaxy in turn. 



3.1.1 Target 1: Jl 03024 ■ 08+ 0524 20.41 

This is the on ly galaxy with spectr oscopic confirmation in 
the sample of IStiavelli et all l)2005l ). We detect a clear (~ 
15cr) emission line (see Fig. [T]) which we identify as Lya at 
ZLy a = 5.973 ± 0.002, in g ood agreement with the redshift 
z-Ly a = 5.970 reported by IStiavelli et all (|2005h . The total 
line flux is Fi^ ya = 22 ± 1.5 x 10 -18 erg s _1 cm" 2 , measured 
by integrating over ~ 23 A (~ 800 km s _1 ) of the spectrum; 
it is n ot possible to compare with the data bv lStiavelli et all 
|2005l ) because these authors did not report a measurement 
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Figure 2. FORS2 spectrum of J103027.98+052459.51. Details as 
in Figure [T] 
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Figure 3. FORS2 spectrum of J103026. 49+052505.14. Details as 
in Figure [l] 



of the line flux. The Lya line flux corresponds to a luminosity 
L-Ly a — 8.6 ± 0.6 x 10 42 erg s _1 in our cosmology, and a 
star formation rate SFR,Lya = 7.8 ± 0.5 Mq yr -1 with the 
conversion (and caveats) given above (see Table [lj. 



3.1.2 Target 2: Jl 03027. 98+052459. 51 

As can be seen from Fig. [21 we detect a weak and nar- 
row emission line which we interpret as Lya emission at 
z = 5.676 ± 0.002. With a line flux F Lya = 4.9 ± 1.1 x 
10 -18 erg s _1 cm -2 , the detection is significant at the ~ 4.5a 
level. Although the line is weak, it does fall in a relatively 
clean region of the spectrum, free from obvious residuals 
from the subtraction of sky emission lines (see top panel 
of Fig. [2j. The corresponding luminosity and star forma- 
tion rate are, respectively, L^ya — 1.7 ± 0.4 x 10 42 erg 
s" 1 and SFR Lya = 1.6 ± 0.4 M yr" 1 (see Table [T]). 
The line has FWHM~ 150 km s _1 after correcting for 
the instrumental resolution. Although this galaxy is rel- 
atively close to the QSO sight-line, at a projected dis- 
tance of 83kpc, no C iv absorption at its redshift wa s re- 
ported by eit her Ryan-Weber. Pettini fc Madaul l|2006l ') nor 
ISimcoe et all (|201lT ). 



3.1.3 Target 3: J103026.49+052505.14 

At a projected distance of 79kpc, this galaxy is the closest 
to the QSO sight-line among the three galaxies targeted by 
our observations. We detect a weak emission line centred at 
Aobs = 8168 A, in a clean region of the spectrum, free from 
strong sky line residuals (see Fig. [2]). We identify the feature 
as Lya emission at z = 5.719 ± 0.002, with a flux -FLya = 
4.3±0.9 x 10~ 18 erg s _1 cnT 2 (~ 4.8cr level detection). This 
corresponds to a luminosity I/Lya = 1-5 ± 0.3 x 10 42 erg s _1 
and star formation rate SFRL ya = 1.4 ± 0.3 Mq yr -1 . 

The emission line redshift differs by Av = —214 km s -1 
from the redshift z abs = 5.7238 ± 0.0001 of the 
strongest Civ absorption system in the spectrum of 



the QSO, with column density Nc 



2.3 x 10 1 



Ryan- We ber. Pettini fc M adau 200J). More recently, 



ISimcoe et al.l 1 20111 ') reported z abs = 5.72438, corresponding 



absorber, measured from near-IR spectra obtained with the 
Folded Port Infrared Echellette (FIRE) instrument on the 
6.5 m Magellan telescope. 

The velocity difference between the Lya emission line 
and the C IV system may be due to random error (~ 2-3<r) 
and/or systematic errors in the wavelength scale of either 
the FORS2 spe ctrum reported here or the near-IR spec- 
tra analysed by Ryan- Weber. Pettini fc Madaul (|2006l ) and 
ISimcoe et al.l (|2011 ). Alternatively, it may be an indication 
that the geometry of the outflowing (or perhaps inflowing) 
gas associated with this galaxy is not spherically symmet- 
ric. The velocity offset could be larger if the Lya emission 
line is redshifted relatively to the systemic redshift of the 
galaxy, as is norm ally the case in galaxies driving larg e 
scale outflows (e.g. iPettini et al.l [200l] ; ISteidel et alj l2010h . 
In any case, as we discuss below, this galaxy is sufficiently 
close to the QSO sight-line for an outflow at the canonical 
«out = 200 km s _1 to have travelled the projected distance 
of 79 kpc in ~ 0.4 Gyr, less than half the age of the Universe 
at z = 5.7. 



3.2 UV luminosities 

In this section we derive the UV continuum fluxes of the 
galaxies from their 2850 magnitudes and use them to deduce 
values of the star formation rate SFRtjv for comparison with 
the values of SFRj^yQ, measured above. At z = 5.7-6.0, the 
ZS50 filter samples the far-UV s pectr um where, according to 
iMadau. Pozzetti fc Dickinson! i|l99Sl) 



L uv (ergs- 1 Hz- 1 ) = 8.0 x lO" 27 SFR uv (M yr- 1 ) 



(2) 



to Av 



-240 km s , and N, 



Civ : 



3.9 x 10 i4 cm"" for this 



for a Salpeter IMF. As mentioned earlier, an IMF with fewer 
low-mass stars, s uch as that appropriate to the Milky Way 
l|Chabrier 1120031 ). would lead to lower values of SFR by a 
factor of ~ 1.8. The conversion in eq. J2} applies to the 
UV continuum at 1500 A, but the intrinsic (i.e. before red- 
dening) UV slope of starburst galaxies is approximately flat 
in F v between 1500 and 1200 A, so that the above conver- 
sion should still apply. On the other hand, the reddening 
correction is unknown for our galaxies; thus the values of 
SFRtjv we derive are strictly lower limits. Furthermore, eq. 
([2} is valid for the ideal case of a continuous star formation 
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Table 1. Lyc? and UV parameters, magnitudes and colours of the targets. 



Object 


Target 1 


Target 2 


Target 3 


R.A. (J2000) 


10 30 24.08 


10 30 27.98 


10 30 26.49 


Dec. (J2000) 


+05 24 20.41 


±05 24 59.51 


±05 25 05.14 


2 850 ( ma §) 


25.94 ±0.05 


26.03 ± 0.06 


25.34 ±0.05 


(«775 - 2850) & 


2.25 


1.23 


1.59 


Z L ya 


5.973 ± 0.002 


5.676 ± 0.002 


5.719 ±0.002 


L h y a (xlO 42 erg s- 1 ) 


8.6 ±0.6 


1.7 ±0.4 


1.5 ±0.3 


SFR Lya (M yr- 1 ) 


7.8 ±0.5 


1.6 ±0.4 


1.4 ±0.3 


F{J V (xl0" 3 ° erg s" 1 cm" 2 Hz" 1 ) 


1.13 ±0.08 


1.29 ±0.08 


2.6 ±0.1 


L uv (xlO 28 erg s" 1 Hz" 1 ) 


6.4 ±0.4 


6.8 ±0.4 


13.6 ±0.7 


SFRuv (M yr- 1 ) 


8.0 ±0.5 


8.5 ±0.5 


17.0 ±0.9 


Physical distance 12 (kpc) 


326 


83 


79 



"Calculated using FLUX.AUTO. 
'Calculated using FLUXJSO. 
c Vacuum heliocentric. 

''Corrected for Lya emission and Lya forest. 
e Transverse distance to the QSO line-of-sight. 



episode lasting for more than 100 Myr; younger ages would 
lead to larger values SFRuv- 

Two corrections that we can apply, knowing the wave- 
lengths of the Lya emission lines and the transmission curve 
of the ACS F850LP filter, are for absorption by the Lya for- 
est and for the contribution to the Z850 magnitudes by the 
Lya emission line itself. For the latter, we simply subtracted 
the flux of the Lya line from the flux measured by SEx- 
tractor in the zsso band (FLUX_AUTO). The correction for 
Lya forest absorption is a multiplicative factor accounting 
for the fact that not all the filter was illuminated in wave- 
length space. We applied the following correction to amend 
for this effect 

rc2 _ rc At 

where /° is the UV flux corrected for Lya emission, At is 
the total area under the throughput curve and Af is the area 
under the same curve for wavelengths A > Aobs(Lya) (thus 
assuming that, to a first approximation, there is negligible 
flux below the Lya emission line). 

As can be seen from Table [TJ SFRuv > SFRL yQ for 
all three galaxies. As mentioned above (Section 3.1), there 
are many possible reasons which can explain the common 
finding that the Lya line luminosity underestimates the true 
star formation rate in galaxies. It is interesting that, out of 
the three galaxies observed, it is the one with the clearest 
evidence for a large-scale outflow (Target 3) that shows the 
lar gest difference betw een SFRuv and SFRL ya - As discussed 
bv lSteidel et al.l (|201ll ). galaxy-scale outflows naturally lead 
to a diffuse, extended halo of Lya emission which is seldom 
captured in its entirety with slit spectroscopy. 

3.3 Color confirmation 

Having confirmed that all three targets are indeed at z = 
5.5-6, we comment bri efly on the colour sel ection for i- 
dropouts. According to Ma lhotra et al. I (120051 ), the colour 
cut 1775 — 2850 > 0.9 in sources with zsso ^ 27.5 will include 
galaxies at z > 5.4 with a completeness of ~ 80%. How- 
ever, galaxies at intermediate redshifts (z ~ 1-2) also have 
1775 — 2850 ~ 1- Therefore, a color cut 2775 — zsso > 1-3 may 



be a stronger criterion for el iminating lower redshift inter- 
lopers, but it was noticed by iMalhotra et all l|2005l) that it 
would lead to higher incompleteness (~ 20-30%) for z ~ 6 
galaxies. 

Concerning our three galaxies, we note that Target 1 
has the reddest colour, 1775 — 2350 = 2.25, but the lowest 
UV flux, once corrections have been applied for Lya forest 
absorption and Lya emission (see Table (TJ. The reason is 
that the Lya emission line makes a significant contribution 
to the flux measured in the zsso filter. 

With zsso = 26.03 ± 0.06 Target 2 is the faintest of the 
three in this band, even though its corrected UV flux is not 
the lowest. Note that with 1775 — zsso = 1-23 this galaxy 
would be missed by a colour selection 1775 — Z850 > 1-3, 
even though it is intrinsically bright: its UV luminosity, 
Luv = 6.8 x 10 28 erg s _1 Hz -1 , corresponds to Luv = 
1.25L* =fi , adopting L * =fi = 5.44 x 10 28 erg s -1 Hz -1 from 
iBouwens et~afl (|2007T l. 

Target 3, which we associate with the C IV absorber at 
z — 5.7238, has the brightest 2:350 magnitude of the three 
objects (zsso = 25.34 ± 0.05 mag) and is also the most lumi- 
nous with Luv = 2.5L* =6 . In this case, the weak Lya line 
does not make a significant contribution to the integrated 
flux in the F850LP filter, resulting in a smaller correction to 
the 2850 magnitude than for Target 1. 

In conclusion, while there is evidence that the colour cut 
2775 — 2850 > 1.3 is effective against low redshift interlopers, 
our (admittedly limited) results show that it can also miss 
a significant proportion of luminous (L > L*) galaxies at 
2 > 5.5. Possibly, a less extreme colour cut, such as 2775 — 
2850 > 1.2, may be more suitable for galaxies at 2 £ 5.7 with 
faint Lya emission. 



4 DISCUSSION 

4.1 Connection between star forming galaxies and 
C iv absorption systems 

To recap, in order to investigate the origin of C iv absorbers 
at 2 ~ 6, we have examined the field of the QSO SDSS 
J1030±0524, whose spectrum shows the strongest C iv ab- 
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sorption doublet in the survey bv lRvan- Weber et all l|2009t) , 
and discovered two Lya emitters at 2 ~ 5.7 close to the QSO 
sight-line (and confirmed a previously known Lya emitter at 
2 ~ 5.97). The redshift of J103026.49+052505.14 (Target 3), 
which is the galaxy closest to the QSO and the most lumi- 
nous of the three, differs by only Av ~ —230 km s~ from 
that of the C IV absorption lineo this finding is highly sug- 
gestive of an association between galaxy and QSO absorber. 

At lower redshifts (z — 2-3), the galaxy-QSO ab- 
sorber c onnection has now been quantified in some de- 



tail (e.g. lAdelben ;er et al.|[20ol . 120051 ; ICrighton et ai]|2010l ; 



ISteidel et al.l 120101 ). In the last of these studies in particu- 
lar, the authors derived a rough relationship between the 
equivalent width of the stronger member of the C IV dou- 
blet, Wo(A1548), and the projected distance between the 
galax y and the QSO sight-line (see Fig. 21 of ISteidel et al.l 
2010). It is thus of interest to verify whether this relation- 
ship also holds at much higher redshifts, at least in the case 
of J103026.49+052505.14. 

The first point to note is that, with a SFRuv = 
17.0 M yr" 1 , J103026.49+052505.14 is typical of the 
galaxies dr iving outflow s at r edshifts z = 2-3 in the 
sample of ISteidel et ail (|2010T ). Second, the rest-frame 
equivalent width of CrvA1548 at 2 a b s = 5.7238 in the 
spectrum of SPSS J103 0+0524 is W (A1548) = 0.65 A 
jRvan-Weber et al l [2009). This value is intermediate be- 
tween the last two points in the Wo(A1548) v s. galactocen- 
tric im pact parameter b plot constructed by ISteidel et al.l 
(2010), between which the line equivalent width drops dra- 
matically from Wo(A1548) = 1.2 ± 0.15 A at b = 63kpc 
to W (A1548) = 0.13 ± 0.05 A at b = 103 kpc. If we 
simply interpolate between these two points, we find that 
Wo(A1548) = 0.65 A is on average expected to be mea- 
sured at impact parameters 6 ~ 70 kpc. For comparison, 
the 13.5 arcsec separation on the sky between the galaxy 
J103026.49+052505.14 and the QSO SDSS J1030+0524 cor- 
responds to an impact parameter b = 79 kpc at z a b s = 
5.7238. Given the uncertainties in the interpolation and the 
likely spread of the average relationship between Wo(A1548) 
and b when applied to any individual galaxy, we conclude 
that the strength of the C IV absorption seen 79 kpc from the 
galaxy J103026.49+052505.14 is typical of the galaxy-scale 
outflows operating at redshifts 2 = 2 — 3. 



4.2 Outflow Speed 

With only the Lya line at our disposal, we do not 
have any information on the speed of the outflow from 
J103026. 49+052505. 14 which has spread metals over radii 
of at least 79 kpc. Furthermore, in the absence of the IR 
photometry necessary to characterise the spectral energy 
distribution of this galaxy, we have no means to estimate 
the age of the starburst. Nevertheless, it is instructive to 

consider some releva nt timescales. 

The analysis bv lSteidel et al.l (|2010T ). together with the 
results of high resolution spectroscopy of g ravitationally 
lensed star-forming galaxies at z = 2-3 (e.g. iPettini et al.l 



3 Averaging the valu es of Av implied by the measuremen ts of 
z aba (C iv) reported bv lRvan- Weber, Pettini fc Madaj l(20Q6tl and 



120021 ; lOuider et all 120091 ; iDessauges-Zavadskv et all |2010| . 

and references therein), has shown that while galactic winds 
may be launched with speed of typically ~ 200 km s~ , the 
gas accelerates outwards reaching speeds of up to 800 km s _1 
before reaching the virial radius. If the bulk of the gas out- 
flowing from J103026. 49+052505. 14 were indeed moving at 
Vout = 800 km s~ , it would reach distances of ~ 80 kpc - 
and could produce absorption in the QSO spectrum - in 
only ~ 100 Myr. 

On the other hand, other authors have considered more 
modest outfl o w vel ocities, largely on theoretical grounds. 
iMartin et al.l |201Q) argued that average outflow speeds 
of less than 200 km s _1 are required to match the line- 
of-sight corr elation function of Civ absor bers at redshifts 
2 = 1.7-4.5. lOppenheimer fc Pavel I 20061 ) favour such ve- 
locities because in their cosmological simulations they are 
found to be sufficiently high to enrich the low density 
IGM a nd yet low enough to avoid overhe ating it. More re- 
cently, iKramer. Haiman fc Madaul (|2010T ) pointed out that 
the rapid rise in the inter galactic density of Civ from 
z = 5.8 and 4.7 discovered bvlRvan- Weber et al.l (|2009l ') and 
I Becker. Rauch fc Sargent] |2009l ') may be inconsistent with 
the near-complete reionization of the IGM by 2 ~ 6 and the 
slower build-up of stellar mass during this epoch. All of these 
observations can be reconciled, however, if there is a delay of 
~ 400-700 Myr between the onset of star formation and the 
appearance of its products (i.e. the metals) at distances of 
~ 100 kpc from the sites of productions. The corresponding 
outflow speeds implied by this argument are also of order 
- 200 km s" 1 . 

Adopting a more modest value v out = 200 km s _1 
for the outflow from J103026.49+052505.14 would imply a 
travel time t ~ 400 Myr to reach the QSO sight-line. In turn, 
this would place the onset of star formation in this galaxy 
at 2 = 8.4 where gala xies are now claimed to be f ound in in- 
creasing numbers fe.g. lBouwens et"al ] |2010l . l201lD . In future, 
deep observations of this and other galaxies in the field of 
the QSO SPSS J1030+0524 at infrared wavelengths should 
allow a rough determination of the age of the galaxy from 
which an estmate of the speed of the outflows may be in- 
ferred. 



4.3 Contribution to the IGM 



ISimcoe et al 



ported I 
. (2011) 



— see Section l3.1.3l 



ISteidel et all (|2010T ) showed that, given the typical extent 
(« 90 kpc) of metal-enriched regions around star-forming 
galaxies at 2 ~ 2, galaxies brighter than ~ 0.3L* can ac- 
count for approximately half of the number of C IV absorp- 
tion systems with VKo(A1548) ^ 0.15 A seen in QSO spec- 
tra at this average redshift. The corollary of this statement 
is that a galaxy with L ^ 0.3L* should be found within 
~ 90 kpc of half the QSO absorbers, with the remainder 
presumably associated with fainter galaxies. 

Given that the properties of the outflow from 
J103026. 49+052505. 14 are at least superficially similar to 
those of star-forming galaxies at 2 = 2-3, as discussed above, 
we can repeat the same exercise at the higher redshifts con- 
sidered here. Specifically, we estimate the number of C IV 
absorption systems per unit of redshift (dN/dz(C iv) pre d) 
expected from galaxies detected in the field of QSO SPSS 
J1030+0524 at 2 ~ 5.7, using the model proposed by 
ISteidel et al] (|2010h . Through a comparison with the ob- 
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served value of dN/dz(C iv), we draw conclusions on the 
properties of galaxies contributing to the metal enrichment 
of the IGM at z ~ 6 and calculate the chances that the 
Zabs = 5.7238 C IV system may be produced by an unde- 
tected, fainter, galaxy instead of J103026.49+052505.14. 

We consider the limiting magnitude in the ACS im- 
ages to estimate an upper limit to the UV luminosity of 
undetect ed galaxies. Ad o pting the z ~ 6 luminosity func- 
tion by iBouwens etafl J20Q3) (M£ v = -20.24 ± 0.19, 
" 3 Mpc" 3 , a = -1.74 ± 0.16) we find that 



= 1.4j£j x 10" 



the ACS images of the field are sufficiently deep to reveal 
galaxies with luminosity L £ 0.9 L% =6 . Integrating the lumi- 
nosity function, we find that the number density of galax- 
ies with L > 0.9L*. =fi is w g a i ~ .287 x 10" 3 h 3 Mpc" 3 . If 
the model by ISteidel et alT l|201Ch is still valid at z ^ 5.5, 
the co-moving cross section for absorption per galaxy is 
Ugai — 1.15 Mpc 2 for an impact parameter b ~ 90 kpc (physi- 
cal). At the 

-^abs — 5.7238 of the C IV absorber, the comoving 
pathlength per unit of redshift is dl/dz ~ 448 Mpc. 

Thus, we want to know the predicted number of absorp- 
tion systems per unit of redshift 

dN dl 

-j— l^lVJpred = n ga l (Jgal — 

with W (A1548) ^ 0.15 A. With the above values 
for rigai, (Tg a i and dl/dz, we find dN/dz(C lv) pro d ~ 
0.15. This estimate is ~ 30% of the observed value 
d NJdzjC iv)o bs = 0.5 ± 0.3 at reds h ift (z) =5.75 estimated 
b v lOppenheimer. Dave Sz Finlatorl ll200Sf). from the results 



of iRvan- Weber et al l (|2009r i. and iBecker. Rauch fc Sargent! 
(2009) for absorption systems with column densities in the 
range logiVrj IV = 14.0-15.0 cm -2 . The shortfall is presum- 
ably higher when w e include weaker Civ systems, since the 
ISteidel et all |2010l ) limit VK (A1548) > 0.15 corresponds to 
column densities logA^Q IV ^ 13.6 cm -2 for unsaturated ab- 
sorption lines, somewhat lower than the column density of 
C IV for which absorption line statistics are currently avail- 
able at these high redshifts. 

We conclude that, if star-forming galaxies at 2 ~ 5.7 
have enriched surrounding regions to radii of ~ 90 kpc, as is 
the case at z = 2-3, it is necessary to reach at least Luv ~ 
0.4 L* =6 (M{j V ~ —19.24, i.e. one magnitude fainter than 
the ACS images considered in this work) in order to detect 
the full population of galaxies linked to absorption systems 
with logJV c ly = 14.0-15.0 cm" 2 . 



5 CONCLUSIONS 

From deep spectroscopy of galaxies in the field of 2 cm = 
6.309 QSO SDSS J1030+0524 we have identified three Lya 
emitters at z ^ 5.7. Of these three galaxies, the brightest 
(Luv = 2.5Z/z=g) and closest to the QSO sight-line (impact 
parameter b = 79 kpc), J103026.49+052505.14, is at a red- 
shift z — 5.719 which differs by only Av ~ —230 km s _1 
from the absorption redshift of a strong C IV system in the 
spectrum of the QSO. 

The chances of a random association between Lya emit- 
ter and C IV absorber are very low. For example, if we first 
consider the idealised case of a random distribution of galax- 
ies with n ga i ~ 1.23 x 10 _3 Mpc -3 , appropriate to galaxies 
of luminosity L ^ 0.4L* =6 , then there is only a 1.1 x 10 -3 



(~ 0.1%) chance of finding such a galaxy within a sphere 
of radius 90 kpc (physical) centred on the QSO sight-line 
at z — 5.7238. However, for this calculation it is more ap- 
propriate to consider a cylinder, rather than a sphere, with 
radius 90 kpc but depth along the line of sight correspond- 
ing to An ~ ±230 km s _1 , the velocity difference between 
the Lya emitter and the C IV absorber. The volume of this 
configuration is approximately six times larger than that of 
a sphere of radius 90 kpc, raising the probability of finding 
a galaxy by chance within this volume to ~ 0.7%. 

A s econd factor to take i nto a ccount is galaxy c lustering. 
Indeed. lStTavelli et all <|2005h and lKim et"afl (|200gh reported 
an overdensity by a factor of ~ 2 at z > 5.5 in the field of 
the QSO J1030+0524. If this average factor applied to the 
environment of the galaxy J103026. 49+052505. 14, the above 
probability would be increased by a factor of two. More im- 
portant, however, is the possibility of significant clustering 
with lower luminosity dwarf galaxies which would remain 
unrecognised given the current limits of the photometry in 
this field. Extending the above calculations further down 
the luminosity function, to galaxies as faint as ~ 0.06 LJ_ 6 , 
would increase the above probability by one order of magni- 
tude. Indeed, the ~ 230 km s _1 velocity difference between 
Lya emission and C IV absorption could be taken as evi- 
dence for an origin of the metals in a dwarf galaxy in the 
same overdensity as J103026.49+052505.14. On the other 
hand, at present we have no way of addressing empirically 
the question of whether such low-luminosity galaxies sup- 
port large-scale outflows of metal-enriched gas, and how 
their mass-outflow rates compare with those of the more 
luminous galaxies in which the outflows are seen directly. 

The result reported here is the first direct observational 
evidence that the circum-galactic medium of early luminous 
star-forming galaxies was already enriched with heavy el- 
ements at 2 ~ 6. Even at these high redshifts, some of 
the properties of the galaxy-wide outflows responsible for 
enriching the circum-galactic medium appear to be super- 
ficially similar to those of their counterparts at z = 2-3, 
which have been characterised in some detail. In partic- 
ular, the strength of the Civ absorption associated with 
J103026.49+052505.14 at a radial distance of 79 kpc fits in 
well with expectations based on data appropriate to lower 
redshifts galaxies, even though the average density of the 
IGM was much higher at z ~ 6 than at z = 2-3. 

We emphasise the importance of extending the avail- 
able photometry of J103026.49+052505.14 to include near- 
IR bands from which the SED of the galaxy can be deter- 
mined with the accuracy required to deduce its age. With 
this information in hand, it should be possible to estimate 
the average speed of the outflow which has transported the 
metals out to at least ~ 80 kpc from the stars that pro- 
duced them. If a speed v out ~ 200 km s^ 1 is appropriate, as 
argued by some, the starburst episode must have lasted at 
least ~ 400 Myr, placing its origin at z ~ 8.4; we now know 
this to have been an epoch when star-formation in galaxies 
was already underway. Higher outflow speeds would imply 
later onsets of the starburst. 

The observations presented here show that it is feasi- 
ble, with dedication, to identify the galaxies linked to QSO 
absorption systems even at the highest redshifts where in- 
tervening metals have been detected in QSO spectra, and 
thus provide an incentive to pursue similar searches in 
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other fields. Finally, we also point out that the colour cut 
1775 — 2sso > 1-3 which has been proposed to avoid low red- 
shift interlopers can miss a significant proportion of bluer 
galaxies at z ^ 5.7 thus affecting mainly the LAE popula- 
tion. We are pursuing this topic in a forthcoming paper. 
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